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}_( ■ The backreaction of nonlinear inhomogeneities to the cosmic expansion is re-analyzed 

' in the framework of general relativity. Apparent discrepancies regarding the effect of the 

^5 \ nonlinear backreaction, which exist among the results of previous works in different gauges, 

. are resolved. By defining the spatially averaged matter energy density as a conserved quantity 

' in the large comoving volume, it is shown that the nonlinear backreaction neither accelerates 

^SJ ' nor decelerates the cosmic expansion in a matter-dominated universe. The present result 

, in the Newtonian gauge is consistent with the previous results obtained in the comoving 

■ synchronous gauge. Although our work does not give a complete proof, it strongly suggests 
^ ' the following no-go theorem: No cosmic acceleration occurs as a result of the nonlinear 

' backreaction via averaging. 

o : 

in ■ 

' The recent observation of the isotropy of the cosmic microwave background ra- 

O , diation (CMBR^ and large galaxy surveys, such as SDSS,^' indicate that the uni- 

\ verse is remarkably isotropic and homogeneous over scales larger than some 100 Mpc. 

' However, it is not straightforward to describe the universe using an isotropic and 

Q^i homogeneous metric, namely, the Friedmann-Lemaitre-Robertson- Walker (FLRW) 

^ \ model, because the local universe is in fact very inhomogeneous. The solution of the 

^ • Einstein equation with an averaged homogeneous matter distribution is not a solu- 

c/3 . tion with a realistic matter distribution, because of the nonlinearity of the Einstein 

\ equation. Thus, it is naturally conjectured that the expansion law of the FLRW 

■ model may be modified by local inhomogeneities. In fact, there have been investiga- 
^ , tions studying this point,'^"^''^''^''"''^''^''^'^''^"'^''"'^^' and some modification has been 

^ \ reported, with apparent discrepancies among these works (for instance, see Refs. ^ 

- - ' and I10|l). The result might depend on the choice of the coordinates as well as the 

definition of the averaging procedure. For these reasons, it has not been possible 
to clearly relate such a nonlinear effect with observations. Furthermore, recent ob- 
servations of Type la supernova^^'^ and the strongly suggest that 
the cosmic expansion is accelerating. Understanding the source of this accelerated 
expansion is one of the greatest unsolved problems in modern cosmology.^'^ This 
acceleration seems to require an unknown type of energy (dark energy, or perhaps 
a cosmological constant). A possible alternative idea to explain the acceleration is 
that the energy resulting inhomogeneities leads to additional terms in the Friedmann 
equation, as if dark energy existed.*^'— Before reaching a definite conclusion re- 
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garding effects due to tlie nonlinear backreaction, we must carefully elucidate the 
construction of a FLRW model with a mean density which obeys the equation of 
state (EOS) for dust matter. Otherwise, we would be misled to conclude that a 
"correction" appears as a deviation from an averaged FLRW model, although the 
averaged matter energy density in such a model does not necessarily satisfy the EOS 
for dust matter. 

Here we attempt to clear up this confusion by carefully constructing an averaged 
FLRW model in the framework of general relativity. Because there is no unique 
choice of the averaged spacetime in the inhomogeneous universe, we take the point 
of view that the averaged density of the matter is conserved in a large comoving 
volume. Our choice seems natural and the most suitable for describing the averaged 
FLRW model with a mean density that satisfies the EOS for dust matter. The 
purpose of this Letter is to show that the nonlinear backreaction neither accelerates 
nor decelerates the cosmic expansion in a matter-dominated universe. It is also 
shown that the backreaction behaves like a small positive curvature term. We show 
that the conclusion does not depend on the definition of the averaging procedure. 

We restrict out study to the case of a vanishing shift vector, i.e., iV* = 0. The 
line element is written as 

ds"^ = -{Ndtf + j.jdx'dx^. (1) 
The unit vector normal to a hypersurface foliated by t = const is denoted by 

71^^= (^,0,0,o) . (2) 

The extrinsic curvature is defined as 

^^•-^7^'7.„ (3) 

where the dot denotes differentiation with respect to time. The Einstein equation is 
reduced to the Hamiltonian constraint, the momentum constraint and (the trace of) 
the evolution equation aJ^ 



(^'i? + {K\) ^ - K'jK{ = WttGE, (4) 

K',i^-K\^=87rGJ.^ (5) 

k\ + NK'jK{ - Tv'j^ = -4ttGN{E + S), (6) 

where '•^^R is the 3-dimensional Ricci scalar curvature, | represents the 3-dimensional 
covariant derivative, and we have 

E = T^,nV = ^Too, (7) 

Ji = -T^^n'' = ^Toi, (8) 

S = Ti,f^. (9) 
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The three-dimensional volume F of a compact domain D on a t = const hyper- 
surface is 



V = [ ^d?x, (10) 



where 

7 = det(7,,). (11) 

Here, V is considered to be a volume sufficiently large that we can assume periodic 
boundary conditions. 

The scale factor a{t) is defined from the volume expansion rate of the uni- 
verse P^'^ 

.« V 1 /■ 1 ., ^ _ ,s. 



Next, we introduce the averaging procedur^ 



With this, we find 



Then we define V^^ as 



{A)^^ I^A^d^x. (13) 



3- = {NK\). (14) 
a 



V) = NK\ - -6%, (15) 



which represents the deviation from uniform Hubble flow. With these definitions, 
one can show {V\) = 0. 

Averaging the Einstein equations, we obtain 



^ = -^{N^{E + S)) + ^{{V\f - V)V\) + \{Nn\ + NK\). (17) 

We assume the matter of the universe to be irrotational dust. Then, the energy- 
momentum tensor is written 

T^"" = pu^'u'', (18) 

where is the four velocity of the fluid flow. 

Up to this point, the treatment is fully general and exact. In order to evaluate 
the R.H.S. of Eqs. (|16jl and l|17|) . we need the actual metric of the inhomogeneous 
universe. We can proceed further either with an exact analytic approach employ- 
ing a class of exact inhomogeneous cosmological solutions or with an approximate 
treatment of the inhomogeneous metric obtained perturbatively. In the comoving 
synchronous gauge, i.e., that with N = l,u^ = (1,0,0,0), an exact analytical treat- 
ment of Eqs. (|16() and H17|) has been carried out by one of the present authors .l^'l^ 
It has been shown that there exists a class of exact inhomogeneous solutions that 
are, nevertheless, homogeneous and isotropic on average, with no backreaction to 
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the cosmic expansion. Approximate approaches in the form of perturbative analy- 
ses have also been performed in several gauges, and some modifications have been 
reported with apparent discrepancies among those works!^'^ In order to clear up 
this confusion, in the following, we re-examine the approximate approach and solve 
Eqs. (|T^ and ((T7|) perturbatively by iteration. 

Let us start from the Einstein-de Sitter background as the zeroth-order solu- 
tion, for simplicity. Next, we obtain the first-order inhomogeneous solution from 
the linearized Einstein equation. It is sufficient to consider the cosmological post- 
Newtonian metric at linear order, '^''^''^'"^ 

= - (1 + 2(l){x)) dt^ + a^{l- 2(l){x)) 6^j dx'dx^, (19) 

where 5ij denotes the Kronecker delta. From the Einstein equation at linear order, 
we obtain 

<P.« = + 2*) , (20) 

^ ^5 = (21) 

where ph is the background density, and contraction has been taken with Sij. Note 
that Eq. (|19|) . together with Eqs. (|2fl|) and (|21j) . represents the first-order solution of 
the Einstein equation in the Newtonian gauge. 

Once we obtain the first-order solution, we can perform the iteration to solve 
Eqs. (|16|) and l|17|) at the next order. Using the first-order solution to the R.H.S. of 
Eqs. ((TB|) and H17() . and retaining up to quadratic order in cp, the averaged Einstein 
equation becomes 

d\ ^ SvrG , 1 , , , , , , 

= ^{Too) + -{^,i<P,^), (22) 

- = -^{Too + Pba^v^) - ;r^(</',^0,.), (23) 
a 3 3a^ 

where v'^ = 6ijV^v^ . 

Equation (|22j) seems to indicate that the nonlinear backreaction expressed by the 
second term on the R.H.S. might increase the expansion rate. However, this is not 
the case. In order to obtain the correct result, we have to determine the mean energy 
density p of dust matter. Most previous works do not give careful consideration of 
this point. The mean density must satisfy 

^ + 3-p = 0. (24) 
a 

Without this consideration, the backreaction cannot be properly related to a devia- 
tion from the Hubble expansion driven by the mean density. Clearly, (Too) does not 
act as p. 

In order to guarantee that p satisfies Eq. pi)) , we have to define the mean density 

as 

P = (Too) + Pba^v^) + ^-i-^(</,^,0,,) = (Too) + (25) 
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where we have used Eq. (|2H) to rewrite the term containing (u^) as that with 
The quantity p is uniquely determined in this form under the following conditions: 
1) it is expressed as a linear combination of the terms appearing on the R.H.S. of 
Eqs. ((TH|l and ((T7I) : 2) it satisfies Eq. pi]). If these conditions are satisfied, then 
Eqs. (|^^ and (^5]) can be rewritten as 

8ttG 1 , , , V , ^ 

-^-7T72('^>*'^-i)' (26) 



aj 3 ^ 9a2 
a AttG 



^ P- (27) 
a 3 

It is thus seen that we can safely call the second term on the R.H.S. of Eq. ()26|) the 
nonlinear backreaction. 

Equations (|26jl and (|27|) reveal the following points. 

• The nonlinear backreaction neither accelerates nor decelerates the cosmic ex- 
pansion. In other words, the cosmic acceleration d/a is determined merely by 
the mean density. 

• The nonlinear backreaction reduces the expansion rate d/a. 

• The nonlinear backreaction is proportional to in the averaged Priedmann 
equation. Hence, it behaves as a (small) positive curvature term in the Pried- 
mann model. This correction might be measured in future space observations, 
such Planck mission. 

It should be noted that Eqs. (PS)) and H27() do not rely on averaging procedures, 
though the explicit form of the mean density does. Por instance, one may simply 
use (A) = f Ad^x^ Changes induced by using a different averaging scheme 
occur simultaneously only inside p of Eqs. (|26|) and H27|) . Purthermore, the nonlinear 
backreaction term is invariant to second order with respect to the choice of the 
averaging procedure. 

It is worthwhile mentioning gauge issues in the backreaction problem. One may 
wonder what happens under a different gauge condition. The comoving synchronous 
(CS) gauge has been employed to investigate the influence of the backreaction on 
quantities used to determine the age of the universJ^ based on the relativistic 
version of a Zeldovich-type approximation.'^''^ Using the CS gauge, one can obtain 
equations similar to Eqs. (|26() and ()27|) . Hence, we reach the conclusion that there 
is no change in the cosmic acceleration, and the averaged Hubble equation has a 
positive spatial curvature as a correction. Therefore, it is strongly suggested that 
the above conclusion is valid for any choice of the gauge. In practice, only the 
post-Newtonian-type and CS gauges have been employed in previous studies of the 
nonlinear backreaction problem, as far as the present authors are aware. With this 
in mind, our conclusion should be considered seriously. 

Pinally, we make a few comments on related works that have appeared recently. 
In §3 of Ref. 22), Ishibashi and Wald give an argument concerning the effect of the 
backreaction, employing an averaging. In spite of their statement that they "point 
out that our universe appears to be described very accurately on all scales by a 
Newtonianly perturbed PLRW metric" , they restrict their analysis to the comoving 
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synchronous gauge. They only consider the "smallness" of possible nonlinear effects 
in the Newtonian gauge. With regard to this point, this Letter gives a new result, 
which shows explicitly that no cosmic acceleration occurs as a result of the nonlin- 
ear backreaction of a Newtonianly perturbed FLRW metric. Hirata and SeljalJ^ 
also reached a negative conclusion concerning the accelerating expansion due to su- 
perhorizon cosmological perturbations. Their argument, however, does not treat 
nonlinear effects due to local (subhorizon) inhomogeneities through averaging, as 
does the analysis carried out in this Letter. 
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